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Abstract. The zero temperature Mott transition region in antifcrromagnetic, spin 
S=5/2 MnO is probed using the correlated band theory LSDA+U method. The first 
transition encountered is an insulator-insulator volume collapse within the rocksalt 
structure that is characterized by an unexpected Hund's rule violating 'spin-flip' 
moment collapse. This spin-flip to S=l/2 takes fullest advantage of the anisotropy 
of the Coulomb repulsion, allowing gain in the kinetic energy (which increases with 
decreasing volume) while retaining a sizable amount of the magnetic exchange energy 
While transition pressures vary with the interaction strength, the spin-flip state is 
robust over a range of interaction strengths and for both Bl and B8 structures. 
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The insulator-metal transition (IMT) in correlated systems (the Mott transition) 
is one of the most actively studied topics in condensed matter systems. [T] The original, 
and most studied, model is that of the single-band Hubbard model (1HM), characterized 
simply by bandwidth W and on-site repulsion strength U. Roughly speaking, for 
U/W > 1 it is an insulator characterized by localized states and local moments, 
while for U/W < 1 the system is a nonmagnetic metal characterized by itinerant 
states. In a more general model, the onset of itineracy would lead to increased 
bonding, hence a volume collapse at the transition. It has recently been emphasized 
that degenerate multiorbital atoms (N orbitals) with multielectron magnetic moments 
behave very differently. The critical interaction/bandwidth becomes (U/W) c ~ y/N 
or even larger, [2] due to the increase in conduction (hopping) channels. Another issue 
is that of a possible orbital selective Mott transition, where only some of the orbitals 
undergo an IMT transition [31 HI \5\ depending on the interactions and the anisotropy of 
the hopping processes. 

A much less studied question, one we address here, is how a multielectron local 
moment disintegrates under reduction of volume. While the moment may be considered 
to be enforced by the strong interaction U (as in the 1HM), the inter-orbital Hund's 
coupling is also a strong factor because it will tend to promote a moment even in 
the itinerant phase. Anisotropic bonding (hopping processes) causes variation in 
bandwidths, and the moment collapse may be orbital selective: some but not all orbitals 
may become spin-paired (doubly occupied), or selected spins may simply flip as the 
kinetic energy overcomes the Hund's coupling but not the Coulomb repulsion. Pressure- 
driven collapses of magnetic signals reported in M 2+ I 2 compounds (M=V, Mn, Fe, Co, 
Ni)[6] and FeO[7], initially interpreted as Mott transitions, have been reinterpreted as 
magnetic (dis)ordering transitions rather than magnetic collapse. [8] In this paper we 
provide predictions for the moment collapse transition in anti-ferromagnetic (AFM) 
MnO at T=0, which proceeds by a different route than any yet envisioned. 

The Mott transition in MnO at room temperature has recently been revealed 
through transport [9] and spectroscopic [10], [H] data at high pressure. Occurring 
entirely within the (spin) disordered phase, there is an insulator-insulator structural 
transformation Bl (rocksalt) — * B8 (NiAs) at 90 GPa, followed by an IMT + moment 
collapse transition at 105 GPa. In fact, this (room temperature) Mott transition at 105 
GPa is the first observed for a 3d monoxide. FeO is reported to remain a magnetic 
insulator to 143 GPa. [8] Unlike for these room temperature experiments where the 
moments are disordered, the magnetic ground state phases we address will be ordered, 
and the resulting symmetry lowering|12] and reduced fluctuations are found to affect 
the character, and probably the mechanism, of the transition. 

Previous theoretical work on MnO at reduced volume has been carried out almost 
entirely in the local spin density approximation (LSDA) and generalized gradient 
approximation (GGA).[13l[Tll[T5l[T2] In these approximations, the small gap at ambient 
volume rapidly closes leading to metallization at much too small a volume. The 
resulting occupation of minority t 2g states at the expense of majority e g states leads 
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to a continuous decrease of the calculated moment well before the volume collapse, or 
moment collapse, transition. Restricted to the Bl structure, GGA gives a metal-to-metal 
moment collapse from 3.5 //# to 1.2 hb at 150 GPa. [13] The main high pressure phase 
is expected to be the B8 (NiAs) structure, [H] and both the crystal and site symmetry 
and structural relaxation have been shown to affect predictions strongly. [15], [12] Fang et 
al. did apply the LSDA+U method to the high volume phase to improve their picture 
of the Mott transition. 

To study the pressure behavior of MnO, we have carried out total energy LSDA+U 
calculations (described below) for both low and high pressure phases. It has recently 
been shown that metal-to-insulator transitions, and even charge disproportion and 
ordering can be modeled realistically with the LSDA+U method. [16] The low pressure 
structure of MnO is well established: it is an antiferromagnetic (AFM) NaCl structure 
with aligned spins in (111) Mn layers, antiparallel with adjacent layers (AFMII). In 
our calculations, we have neglected the small rhombohedral distortion angle. There are 
two simple arrangements of the Mn spins in the B8 phase, ferromagnetic (FM) or AFM. 
Calculations, including structural optimization, show the AFM phase to be energetically 
favorable by 0.2 eV per formula unit, for a wide range of pressures. 

Results we present below use the LSDA+U methodjTT] in the rotationally invariant 
form[18] as implemented in version 5.20 of the full-potential local orbital band structure 
method (FPLOpU [20]). Q). The implementation of the LSDA+U method in this code 
has been provided in detail by Eschrig et a/.[2l] The all-electron aspect of this code is 
important, since even small-core pseudopotentials cannot reproduce all-electron results 
under volume reduction. [25] The full-potential aspect can be important also, on the 
oxygen site as well as on the Mn site, and the non-spherical aspect of the potential will 
grow as the volume is reduced. Due to the unexpected nature of the reduced spin state, 
several results were checked, and reproduced, using the Wien2k code.[26j LSDA+U gives 
two distinct spin states for both Bl and B8 structures: a low pressure high-spin (HS) 
S = | configuration and a high pressure, low spin (LS) S = \ state. The equation of 
state (EOS) curves are displayed in FigJTJ From the enthalpies we obtain a first-order 
magnetic transition from HS-B1 to LS-B1 at P cl =123 GPa, followed by a structural 
transition to LS-B8 at P C 2=130 GPa. The isostructural volume collapse at P c \ is « 5%. 
However, these results vary with the choice of U and J (we use 5.5 eV and 1.0 eV 
respectively), which is discussed later. This part of our results has been compared 
with those of other correlated band theory results recently. [27J An unusual feature of 
the present results is the persistence of the LSDA+U bandgap up to higher pressures, 
beyond the observed IMT at room temperature. [9] Thus the magnetic and structural 

\ A single numerical basis set for the core states (Mn Is2s2p and O Is) and a double numerical basis 
set for the valence sector including two 4s and 3d radial functions, and one 4p radial function, for Mn, 
and two 2s and 2p radial functions, and one 3d radial function, for O was used. The semi core states 
(Mn 3s3p) are treated as valence states with a single numerical radial function per nZ-shell. We have 
used the strong local moment form of the LSD A double-counting correction fT?] [211 E21 [23] that has 
been become known as the 'atomic-limit' form. [23] The Slater parameters were chosen according to 
U = F = 5.5 eV, J =±(F 2 + F i ) = l eV and F 2 /F A = 8/5. 
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Figure 1. The calculated total energy/MnO versus volume for the low pressure (NaCl) 
and high pressure (NiAs) structures of MnO. The filled symbols denote the calculated 
energies and the continuous lines are the least square fitted curves to the Murnaghan 
equation of state for high and low spin configurations respectively. The inset clearly 
elucidates the order of the transitions, NaCl (high-spin) — * NaCl (low-spin) — * NiAs 
(low-spin). 
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Figure 2. LSDA+U DOS, at the indicated volumes, projected onto symmetrized Mn 
3d orbitals in (left panels) the rhombohcdral Bl AFMII phase and (right panels) the 
B8 AFM structure. In each case, the top subpanel is for the high spin state, while the 
bottom is for the low spin state. See the text for the definition of the a g ,eg,e^ labels. 
The overriding feature is the spin-reversal of the m = ±1 e g orbitals between the two 
volumes. Broadening of the a g states in the LS-B8 DOS is due to the direct Mn-Mn 
d-overlap in the z direction in B8 structure (Mn lies on a simple hexagonal sublattice) . 



transitions we discuss are always insulator-to- insulator, and for the magnetically ordered 
state at T=0. 

To help in understanding the mechanism of the transition, the densities of states 
(DOS) projected onto each of the 1 = 2 irreducible representations are displayed in Fig. 
[21 referenced to the rhombohedral axis of the Bl AFMII phase, and equivalently the 
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hexagonal axis of the B8 structure. Not evident from this figure are two potentially 
important features. [28] First, the 'charge transfer' energy increases as the volume 
decreases; that is, the O 2p levels drop in energy relative to the Mn 3d states, reflecting 
an increased tendency toward the fully ionic limit that competes with the increased 
hybridization as the Mn and O ions approach each other. Secondly, the crystal field 
splitting between E g and T 2g states increases under pressure, which competes with 
Hund's exchange and with correlation effects. 

In each case the Mn site symmetry splits the 5 d-orbitals into two doublets e^ m ', 
|m| = l, 2, with e x g — > {xz, yz}; e 2 g — > {x 2 — y 2 , xy}, and a singlet a g — > 3z 2 — r 2 (m = 0). 
The actual state realized at a given volume or structure is characterized by two e g pairs, 
obtained by unitary mixing, schematically: 
e a g = cos/3 e) -sin/5 e 2 g , 
e b g = sin (3 e x g + cos (3 e 2 . 
We have observed and quantified the mixing angle (3 versus pressure, but near the critical 
pressure it simplifies to /3 ~ 0. 

The HS states in both Bl and B8 structures are simple - each 3d orbital is filled once 
with spins aligned leading to an S=5/2 spherical ion. The LS DOS for both structures 
reveal the essence of the HS-LS transition: the LS state is obtained by simply flipping 
the spins of the e x g orbitals. This result shows how the LSDA+U method differs in an 
essential way from LSDA, where the moment decreases continually with volume. [13] 
metallization occurs at low pressure, and decrease of the moment implies rapid collapse 
of the exchange splitting, resulting in doubly occupied orbitals with zero net spin. In 
this LS state, each of the 3c? orbitals remains singly occupied, the charge density remains 
spherical while the spin density becomes highly anisotropic, as illustrated vividly in Fig. 

El 

To understand the origin of this spin-flip state we have performed an analysis of 
the LSDA+U method. The flavor we have used is the "atomic limit" (AL).[2Tl l22l [23] 
We consider first the pressure-induced change in kinetic energy relative to the potential 
energy of the HS and LS states. We note first that the states, whose spins flip in 
the LS state, are \E g and only \T 2g in terms of the cubic states we are more familiar 
with. For the e 2 pair this ratio is opposite, and close to the overall mean (f-E^, \T 2g ). 
It is the cubic E g states (and consequently the e g pair) that have the strongest (dpa) 
overlap with O ions, and thus are most affected by pressure and give the greatest gain 
in kinetic energy. 

Now we consider the effects of the very large anisotropy of the LS Mn ion. To 
separate the effects of U from those of J, we split the energy expression into the isotropic 
interaction, and the remaining anisotropic part: E AL = E amso + E tso . The isotropic part 
reads 

E lso = - (U - J) i nS (1 - nS )) > °, (2) 

where n s is the spin-dependent occupation number matrix of the 3c? shell, and Tr denotes 
a trace of the orbital indices. It is easily seen that E lso takes its minima for integer 
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Figure 3. Isosurface plots of Mn ion spin density, with red and blue shading indicating 
opposite sign. Left: before collapse, showing the spherical S=5/2 ion, with (111) 
layers of aligned spins. Right: after collapse, revealing the anisotropic S=l/2 ion. The 
magnetic order remains AFMII. 

occupations, in which case we get E lso = 0. However, the resulting potential matrix, 
to be added to the Kohn-Sham Hamiltonian, is non-zero and has the effect of lowering 
occupied orbitals by — \ (U — J), while raising unoccupied orbitals by ~ (U — J). This 
separation stabilizes insulating magnetic solutions, and also favors HS exchange energy 
contributions from the magnetic part of the LSDA functional. This tendency opposes 
the observed HS-LS transition, but operates equally independent of volume. The action 
of the isotropic term is often the dominating effect of the LSDA+U method, but the 
common practice of discussing the effects of the LSDA+U method in the isotropic limit 
misses the physics of this transition, as we now illustrate. 

The anisotropic term reads (in the representation in which n s is diagonal, for 
simplicity) 

E amso = lj2J2< [A£W - S ss ,AJ m ,] <. (3) 
ss' l^rn 

The interaction matrix elements are defined as AU m)1 = w™ 1 ^ — U, AJ mfl = — J — 
(U — J) 5 mfJi in terms of the full matrix interaction w™*t, and these differences do not 
contain U (whose effect is included entirely in the isotropic term). These differences 
describe pure anisotropy, summation of either index of either one gives a vanishing 
result. Thus a filled spin subshell will not contribute to E amso , as expected intuitively. 

The HS state is favored by the LSDA spin polarization. However, the LS state 
still has all fully polarized orbitals, so the energy difference will be much smaller than 
for usual LSDA S=5/2 and S=l/2 moments. The resulting ratio of exchange energies 
can be estimated from the integral over the square of the spin density, which gives a 
value Ei/ 2 /E 5 / 2 ~ 0.31. This reduction of exchange energy in the LS state is much less 
dramatic than the estimate from the simple formula E x = — Jst ° nci M 2 . A more detailed 
explanation of this energy ratio is given in the Appendix. 
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The LS configuration requires the a g orbital to be singly occupied. The remaining 
4 electrons then are distributed in pairs among the e a g ,b doublets. Analysis shows that 
the two occupation patterns for which the electrons doubly occupy either the e g or e 2 g 
orbitals have the same anisotropy energy of about E aniso > — 0.3J; a spread of energies 
arises from allowed mixing of e g symmetries. The resulting spin density of such states 
has a 9 -derived shape. 

The remaining two patterns are obtained by occupying e a g in the up-channel and 
e b g in the down-channel or vice versa, both of which results in strongly anisotropic spin 
densities. The lowest energy of E amso m — 1.85J is obtained when the e 2 g and a g are 
occupied in the same spin channel, while e g is occupied in the opposite channel (note, the 
mixing angle (3=0). The other solution (e l J || aj) has E aniso « —1.14 J. The dependence 
of the energies on the mixing angle complicates the discussion. 

However, there is a gap of ~ 0.84J between the spin-flipped and non spin-flipped 
solutions, which is not closed by any mixing. It turns out that this gap is due to the 
density-density anisotropy (AU in Eq. (J3])), which pushes the (density-) non-spherical 
non spin-flipped patterns up in energy, while it is zero for the (density-) spherical 
spin-flipped patterns. The exchange and self-interaction-correction contributions to the 
anisotropy (A J in Eq. (151)) is nearly of the same size for the non spin-flipped and 
spin-flipped occupation patterns, and hence not changing the energy separation of these 
pattern classes. However, it further discriminates the two spin-flipped patterns. 

The anisotropy (which is solely controlled by, and proportional to, J in the 
LSDA+U method) of the interaction favors an occupation pattern which maximizes 
the spatial distance between the electrons (under the constraint of S = ~), while the 
isotropic term oc U — J merely selects insulating over metallic solutions. We indeed find 
that both the spin-flipped and non-spin flipped solutions may be found in LSDA+U 
calculations, separated by an energy of the order derived here. When J is decreased 
to zero, the energy difference between these solutions shrinks, leaving only the LSDA 
anisotropy energy difference, which turns out to be very small: about 0.3 eV at ambient 
pressure and decreasing to zero close to the transition. Of course the LSDA part of the 
functional includes some anisotropy effects, however the main difference in the influence 
of the LSDA-anisotropy and the J-anisotropy Eq. ([3]) is its action on the Kohn-Sham 
states. The LSDA potential contributions act on all states, while the LSDA+U potential 
matrix acts orbital selective. It is this very selectivity, which makes the whole LSDA+U 
machinery work, by mimicking the suppression of occupation number fluctuations due 
to correlations. In the same way as the Hubbard band split is not attainable in LSDA, 
the anisotropy effects are largely suppressed. This suppression is nicely confirmed by the 
observation of the vanishing energy difference between the two LS solutions as described 
above. 

With increasing pressure the kinetic energy gain becomes more and more 
competitive with the exchange energy due more to the increasing crystal field splitting 
than to the bandwidth. This competition usually leads to a (partial) collapse of the 
magnetic moment. LSDA calculations give a moment of about 1.5 \xb at our transition 
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pressure, which is far from the HS value of 5 /jLb thus clearly showing that at the 
transition pressure first Hund's rule is strongly suppressed. In the LSDA+U method, 
the isotropic term forces the HS solution to have full spin-moment, while the LS solution 
allows a larger gain in kinetic energy, hence bringing along a transition from HS to LS 
at some pressure. The anisotropy contribution of LSDA+U is zero for the HS state 
and negative for the flipped LS state. This will further lower the energy of the flipped 
LS state against the HS state, resulting in a lower transition pressure. Moreover, this 
anisotropy contribution is smaller for a non spin-flipped solution, which rules out this 
solution. Since the anisotropy term offers a way to keep a sizable amount of magnetic 
exchange energy, while gaining kinetic energy, it is this "unusual" state, which is realized 
after the transition. 

We have described here how LSDA+U energies for MnO under pressure predict 
an unexpected mode of collapse of the Mn moment at zero temperature: each of the 
3d orbitals remains polarized, and an S=5/2 to S=l/2 reduction arises from a simple 
spin-flip of the symmetry-determined e l g doublet that has the strongest overlap with 
neighboring O 2p orbitals. This S=l/2 moment in the high pressure phase is consistent 
with the interpretation of xray emission data by Rueff et a/.[H]; Yoo et al. were less 
specific about the value of the (clearly small) high pressure moment but presumed total 
collapse. The partial spin-flip collapse obtained here occurs in both the Bl (rocksalt) and 
B8 (NiAs) structures, calculated by two different codes, and occurs at similar volumes. 

The transition we find is first-order and insulator-to-insulator, both of which 
insinuate the smallness of fluctuation effects and make the LSDA+U approach an 
appropriate one. Due to the imprecisely known values of U and J there is an associated 
uncertainty for the calculated equation of state, which depends on the values chosen. 
The functional dependence of these interaction energies (U, J) on the density is not 
known and we have neglected the volume dependence. The Wien2k code [26] includes 
a constrained LSDA algorithm that enables calculation of U, resulting in a value of 
6.5 eV for both for the equilibrium volume and close to the transition. To identify 
the range of variation, we performed EOS calculations for various values of U and 
J. The resulting transition pressures are depicted in Fig. HI As should be expected 
from the importance of the correlation corrections, the transition pressure P C (U, J) is 
quite dependent on the parameters. For all values of U, P c decreases with increasing 
J. The strong variation confirms our analysis above, which identified the anisotropy 
(proportional to J) as crucial in determining the ground state at a given volume. The 
energy separation of the HS and LS curves decreases with J, and this change decreases 
the transition pressure (the shape dependence of the energy curves on J is minor). 
On the other hand, the potential matrix element effects (which shift the corresponding 
eigenvalues) tend to increase the Hubbard splitting with increasing U, which moves to 
stabilize the HS solution against the LS solution and leads to a monotonic increase of 
P c with U. Note that for reasonable values of J— 0.6-1 eV, the dependence on U lessens 
compared to stronger variation for unreasonably small values of J. The very much to 
large transition pressure for J = eV compared to experiment is a strong argument that 
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Figure 4. Dependence of the high-spin to 'spin-flip' low-spin transition on the 
parameters, U and J. 

the anisotropy effects are not sufficiently described by the LSDA part of the functional 
and hence have to be included in an orbital selective manner. 

Although the occupation number fluctuations at this transition should not be a big 
factor, as the volume is reduced and the bandwidth increases, these fluctuations will 
tend to increase. The effect can be modeled by adopting a smaller U than the ambient 
pressure value; hence the value of 5.5 eV that was used for most results presented here 
becomes justified. On the other hand, the value of J is only weakly screened by the 
environment, and it is sensible to chose a volume- independent value that resembles 
the atomic/ionic situation, where one usually finds J =0.7-1 eV for transition metals. 
Altogether, the calculated transition pressures vary between 100 and 170 GPa for a 
reasonable choice of J, which is quite satisfactory given the uncertainties of the LSDA+U 
approach. One should also keep in mind that the zero temperature transition is expected 
to occur at a higher pressure than at room temperature. 

However, we want to stress our main point. Although the LSDA+U method is not 
capable of precise predictions of P c due to the uncertainties just discussed, the spin- flip 
character of the predicted ground state is highly stable against changes of the interaction 
parameters. For all parameter sets corresponding to the data points in Fig. H]we obtain 
the spin-flip LS solution as the ground state. 

Now we summarize. For MnO at T=0 in the AFM ordered phase, an unusual 
moment collapse S = 5/2 — ► 1/2 is predicted before the Mott transition (metallization, 
or itineracy of the 3d states). These results are robust: the spin-flip state is obtained 
by two different codes, and for a substantial range of choices of U and J. The value 
of U = 5.5 eV used here corresponds to U/W ~1.5 in terms of the full 3d bandwidth 
W. The role of J is central to this transition, but in an unexpected way. Hund's first 
rule, which is encouraged by the spin-exchange aspect of J, is violated at the transition, 
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whereas the anisotropic Coulomb repulsion that is proportional to J becomes the driving 
force. Together with an orbitally-dependent increase in kinetic energy, the result is an 
orbitally-selective spin-flip collapse of the moment at an insulator-to-insulator transition. 

The order and type of transitions under pressure we obtain differ from that 
observed at room temperature. [HI [10] It is established however that structural phase 
boundaries can be strongly temperature dependent in transition metal oxides, [29] so 
there is no contradiction. The predicted pressure range is accessible to diamond 
anvil cell experiments, and the ordered-phase moment can be probed by Mossbauer 
spectroscopy. [6] 
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Appendix 

The modern flavour of LSDA+U explicitly excludes first Hund's rule from the expression 
added to the functional, arguing that this contribution is dealt with better within the 
LSDA part of the functional. In order to obtain estimates of this LSDA-contribution, 
we will derive an approximate expression for the LSDA xc-energy E xc suitable for the 
situation we discussed in this paper. 

We expand E xc up to second order in variations of the density. A natural choice for 
a reference density would be a spherically averaged non-spin-polarized density around 
the atom center. We denote the reference density by p . Our interest is in the effect of 
different orbital occupations on the magnetic ion. We can describe the spin-density of 
the /-shell by 

P S W=E^( r X,nC( r ) > 

mn 

where n s mn denotes the generalized occupation number matrix for spin s — ±1, 0/ m (r) 
are suitable orbitals and the indices m, n run over the orbitals of the shell. In the usual 
manner one can write the orbitals as fixed atom-like functions 

<pim(r) = Ri(r)Yi m (r) , 

putting the flexibility into the occupation number matrix. We can introduce the 
occupation number matrix n = J2 S U S f° r the charge density p(r) = J2sP s ( r ) an d the 
occupation number matrix m = J2 S s lk s for the magnetization density m(r) = J2 S sp s (r). 
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The particle number of the spin channel s is N s = Trn s , which gives the total number 
of particles as N = J2s N s and the magnetic moment as M = Y^ s sN s . The spherical 
and spin average of p s is included in the definition of the reference density, hence the 
density variation due to different occupation patterns is 

5p s (r) = [R(r)] 2 J2Y lm (r)Sn^ n Y; n (f) , (A.l) 

ran 

where the variation of the occupation numbers <5n s is measured with respect to the 
averaged occupation numbers (A = 21 + 1) 

n mn,0 ~ ®rnn ■ 

Since the reference density is non-polarized the variation of the magnetization density 
equals the magnetization density itself: m = 5m, m(r) = Sm(r). 

To keep things simple, we restrict the discussion to the local-density form of the 
xc-energy 

E xc = J pe(p + ,p~)d 3 r . 
The second order variational expansion around the reference density then reads 
E xc = E XC)0 + J V xcfi {r)5p{r)d 3 r + J B xcfi {r)m{r)d 3 r 

~P (r)(5p(r)) 2 + 2Q (r)(5p(r)m(r)) + K (r)(m(r)f 



d 6 r . 



1 

+ 5 

The xc-potential and the second order xc-kernels are spherical due to our spherical 
reference density (T4c,o( r ) = Vxc,o(r)) and the xc-field is zero, since the reference density 
is non-polarized. Using this information and the shape of the density variation Eq. 
(lA.lj) we arrive at 



E xc = E xcfi + v 5N 



2 

mnm'n' 



with the variation of the particle number SN = Tr 5n, with the angular coefficients 

a m n, m 'n' =4:71 Yi m (f ) Y\ n (f ) Y\ m i (f ) Y[ n i (r)dQ (A. 2) 



and with the radial integrals 

v = f ^ C)0 (r)[i^(r)]Vdr 



k = ±- I K {r)[Ri(r)] 4 r 2 dr (A.3) 



4.TT 

Po = ^ I P Q (r)[Ri(r)] 4 r 2 dr 



% = ^I Qo( r )Wr)] 4 r 2 dr . (A.5) 



For I < 2 and real spherical harmonics we obtain 

A 



"mn,m' n' 



[8mn8m'n' + 25 mn i6 nm >] , (A-6) 



A + 2 
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Po(5Ny + 2q Q 5NM + k M 2 



which leads to the simple expression 

E xc = E xcfi + v SN 
A 

+ 2(A + 2) 

2p Tr(5n) 2 + 4g Trunin) + 2k Tr(m) 2 ; 

In discussing the MnO case, we consider only 3d 5 occupation patterns, hence 5N = 0. 
For a spherical magnetic occupation pattern we have to set 5n s mn = 5 mn ^-, 5n mn = 0, 
8m mn = 8mn%, resulting in SE XC = ^M 2 , which suggest the interpretation k = — | 
with the Stoner parameter I. 

For the HS and spin-flipped LS pattern (SF-LS) we get <5n = 0, since these patterns 
correspond to a spherical charge density. The magnetic occupation numbers are diagonal 
and equal to ±1, hence Trm 2 = A. So we get 

5E™ = --M 2 = --25 

J - L 4 4 

,^sf ls I nr2 I 2 2 I J48 J55 

= — M 2 ( A Tr m 2 - M 2 ) = — 1 = . 

4 4A+2V — J 4 47 47 



The SF-LS energy contains a large contribution, which is related to the non-sphericity 
of the spin-density. It accounts for pa 7/8-th of the whole xc-energy of this configuration. 
The ratio between the energies of these two configurations is E SF _ LS / En S = 11/35 ~ 
0.31. 

For the non spin-flipped pattern (NSF-LS) the diagonal occupation number 
matrices read Sn = (1, —1, 0, —1, 1) and 8m = (0, 0, 1, 0, 0), which gives for the change 
in xc-energy 

-.NSF-LS _ ^71^2 I ^ / , A m_/_\2 » ^ Po 2A 



^nsf-ls = _ i M 2 _ ATr(m) 2 - M ) - -™-^Tr(5n) 2 

4 4 A + 2 V v — 7 / 4 k A + 2 v ~ ; 

/ J8_/p 40 
~ ~ 4 ~ 4 7 ~ 4A^Y ' 
Again, the second term is due to the non-sphericity of the spin density. However, it is 
much smaller than for the flipped case. The third term, proportional to po is related 
to the non-sphericity of the charge density of the shell. Estimates of po from actual 
calculations give a value of po ~ |/co, hence the third term is roughly of the same size as 
the sum of the first two terms: SE^ F ~ LS pa — 15, which is only two-thirds of the energy 
of the spin flipped case. 
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